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I.  Introduction 

\ 

-'Since  the  solution  of  the  KdV  equation  by  Gardner,  Greene, 

Kruskal  and  Miura*  in  1967,  the  inverse  scattering  method  which 

they  discovered  has  been  extended  to  solve  a  large  class  of  nonlinear 

PDE's.  Almost  all  such  solutions  have  been  discovered  either  by 

inspired  guesswork  or  by  starting  with  a  scattering  problem  and 

working  backwards  to  see  what  PDE's  it  could  be  used  to  solve.  However, 

a  method  for  seeking  an  inverse  scattering  solution  to  a  given  equation 

2 

has  emerged  from  the  work  of  Wahlquist  and  Estabrook. -  They  have 
studied  integrable  nonlinear  PDE's  in  terms  of  differential  forms. 

Their  method  has  been  reformulated  in  purely  classical  terms  by 
Corones  and  further  developed  as  a  useful  computational  tool  by 

34  5 

Corones  ,  Kaup  ,  and  Newell  . 

The  Wahlquist-Estabrook  (WE)  method  systematically  associates  a 
set  of  Lie  brackets  with  a  given  nonlinear  PDE.  If  this  Lie  algebra 
can  be  closed,  consistent  with  the  Jacobi  identity,  then  an  associated 
linear  representation  gives  the  inverse  scattering  transform,  when  one 
exists.  Even  when  an  inverse  scattering  transform  does  not  exist, 
if  the  Lie  algebra  can  be  closed  the  associated  linear  representation 
can  provide  useful  information  about  the  solutions  of  the  nonlinear 
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PDE.  Indeed,  Kaup  has  shown  that  a  WE  analysis  of  Burger's  equation 
leads  to  the  solution  by  a  Cole-Hopf  transformation. 

By  far  the  most  laborious  task  in  applying  the  WE  method  to  a 
complicated  nonlinear  PDE  is  that  of  closing  the  Lie  algebra. 
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^'^Ihis  task  Involves  repeated  application  of  the  Jacobi  identities,  which 
in  turn  requires  the  evaluation  of  a  large  nurrber  of  Lie  brackets  which 
must  be  looked  up  in  a  table.  The  simple  but  repetitive  nature  of  this 
task  suggests  that  we  attempt  to  mechanize  it  using  a  computer.  However, 
we  must  also  require  that  our  program  automatically  perform  algebraic 
simplifications  at  each  step:  properly  distributing  the  Lie  bracket  over 
a  sum  of  terms,  combining  like  terms,  multiplying  numerical  coefficients, 
etc.  We  have  therefore  used  the  symbolic  manipulation  language  MACSYMA 
to  implement  our  algorithms.  The  language  has  built-in  capabilities  for 
performing  the  standard  simple  algebraic  manipulations  from  which  we 
build  up  our  algorithms.  We  also  make  use  of  MACSYMA* s  built-in  program 
for  solving  simultaneous  ^Xnear  equations  algebraically. 

Closing  the  Lie  algebra  ordinarily  involves  sane  amount  of  guesswork. 
To  effectively  use  the  conputer  we  must  first  develop  a  well-defined 
algorithm  which  mechanizes  as  much  of  the  task  as  possible.  We  are  willing 
to  pay  for  this  increase  in  mechanization  by  increasing  the  amount  of 
algebra  that  must  be  performed.  However,  we  also  require  that  our  imple¬ 
mentation  be  sufficiently  flexible  to  allow  us  to  make  some  use  of  the 
insight  we  get  as  our  algorithm  progresses.  We  describe,  in  this  paper, 
our  development  and  implementation  of  an  algorithm  which  serves  these  pur¬ 
poses.  We  also  present  an  exanple  of  its  use  to  find  previously  unknown 
solutions  to  a  phsyically  interesting  set  of  coupled  PEE*s. 

The  PDE's  which  we  take  for  our  example  are  known  in  plasma 
physics  as  the  wave  kinetic  equations  Cor  weak  turbulence  three-wave 
equations) : 
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They  describe  the  interaction  of  three  wave  packets  which  are 
slowly  varying  on  the  space  and  time  scales  of  the  corresponding 
wavelengths  and  periods,  yet  sufficiently  broad  in  Fourier  space 
to  justify  the  use  of  a  random  phase  approximation.  In  Eqs.  (1) 
the  N' s  are  the  mean  squared  amplitudes  of  the  waves  and  the  v^ ' s  are 
the  group  velocities.  The  solutions  to  Eqs.  (1)  ignoring  spatial 
derivatives  (time-only)  are  given  in  standard  plasma  physics  texts6. 
The  time-only  equations  are  widely  used  to  describe  wave-wave 
interactions  in  a  turbulent  plasma.  A  more  complete  description 
of  our  results  concerning  Eqs.  (1)  is  given  elsewhere  . 
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II.  The  Wahlquist-Estabrook  (WE)  Method 

We  give  a  brief  outline  of  the  WE  method  here,  sufficient  for 
explaining  the  role  of  Macsyma  in  our  work.  The  workings  of  the 
method  should  become  somewhat  clearer  when  we  go  through  the  example 
in  Section  V.  A  more  detailed  description  of  the  WE  method  is  given 
in  the  references  cited  in  the  introduction. 

The  basic  idea  of  the  method  is  to  imbed  a  given  nonlinear  PDE, 

^  *  F (u,Uj^,u^,  • . . ) ,  (2) 

as  the  consistency  condition  for  a  pair  of  linear  equations: 


p^V'w***)  t* 

(3a) 

Q(u,ux,...)  <fr. 

(3b) 

Here  <f>  is  an  auxiliary  vector  function  of  x  and  t,  and  P  and  Q  are 
matrix  functions  of  u  and  a  finite  number  of  its  derivatives.  The 
condition  for  consistency  of  Eqs.  (3a)  and  (3b)  is 
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or 
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(5) 
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If  P  and  Q  are  given  functions  of  u,  (5)  is  a  PDF.  (in  general  nonlinear) 
for  u.  We  would  like  to  choose  the  matrix  functions  P  and  Q  so  that 
Eq.  (5)  is  the  same  as  Eq.  (2). 

Applying  the  chain  rule  to  Eq.  (5),  using  Eq.  (2)  to  replace  ut 
by  F,  and  defining 

[P  »Q]  s  PQ  *  QP  (6) 


we  obtain 


PuF 


\Uxt 


-  Vx  -  Vu: 


XX 


+  [P»Q]  -  o. 


(7) 


Note  that  the  highest  derivative  of  u  appearing  explicitly  in  (7)  is 
not  an  argument  of  P  or  Q.  Since  the  initial  value  of  u  is  arbitrary, 
we  therefore  demand  that  the  coefficient  of  this  term  vanish.  Solving 
the  resulting  equation,  we  substitute  back  into  (7),  and  then  we 
iterate  this  process.  It  will  be  clearer  how  this  works  when  the 
example  is  presented  in  section  V.  When  this  procedure  can  be 
completed  we  are  led  to  a  solution  for  P  and  Q  of  the  form 


l 


p*  l  X.(u,u  ,...), 
j-1  1  x 

rt 

oo 

w 

m 

Q  *  L  Y.(u,u  ,•••), 
j-i  3  x 

OO 

cr 

w 

r 


■ t 


where  various  of  the  commutators  between  X^'s  and  Yj's  are  determined. 

If  we  can  close  the  "unclosed  Lie  algebra"  determined  by  the 
Xj ' s,  the  Yj's,  and  the  commutation  relations  we  have  found,  then  we 
can  find  a  linear  representation  for  the  Lie  algebra.  This  linear 
representation  then  gives  us  a  solution  for  the  matrices  P  and  Q. 

For  those  nonlinear  equations  which  are  soluble  by  means  of  an  inverse 
scattering  method,  a  free  parameter  appears  in  this  closure.  This 
free  parameter  corresponds  to  the  eigenvalue.  For  nonlinear  equations 
soluble  by  a  Cole-Hopf  type  of  transformation,  no  such  free  parameter 
appears . 

Closing  the  Lie  algebra  is  a  critical  step,  often  involving 
much  guesswork  and  much  laborious  algebra.  Systematic  procedures 
that  work  for  relatively  simple  equations  like  the  KdV  equation  become 
discouragingly  laborious  when  applied  to  the  Lie  algebras  for  more 
complicated  equations.  We  therefore  use  Macsyma  to  perform  the 
algebra  involved  in  this  step.  Two  algorithms  which  we  have  implemented 
on  Macsyma  for  this  purpose  are  described  in  the  next  section.  Taken 
together,  the  two  algorithms  appear  to  provide  a  very  general  method 
of  closing  the  algebra. 

We  should  mention  here  the  implementation  by  Wahlquist  and  Estabrook 
of  the  method  of  differential  forms  on  MACSYMA. 8  While  their  concern  has 
been  to  bring  to  bear  the  full  power  of  the  formalism,  we  use  only  the 
properties  of  the  Lie  bracket.  We  have  gone  beyond  their  work  in 
developing  and  implementing  the  set  of  algorithms  described  in  this  paper 
for  closing  an  "incorrplete"  Lie  algebra. 
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III.  Algorithms  for  Closing  the  Algebra 

Systematic  procedures  for  closing  the  Lie  algebra  involve  at  each 
step  the  evaluation  of  the  Jacobi  identity 

[A,  [B,C]  ]  +  [B,  [C,A] ]  +  [C,  [A,B]]  =  0  (9) 

for  all  triplets  A,B,C  such  that  the  commutator  is  known  for  each 
pair  of  the  triplet.  The  identities  obtained  may  determine  previously 
unknown  brackets.  If  they  do,  then  the  new  brackets  are  substituted 
back  into  the  Jacobi  identities.  The  process  continues  until  no  more 
new  commutators  are  determined.  Mot  only  is  this  an  extremely  tedious 
process  if  there  are  a  large  number  of  elements,  but  it  is  also  a 
process  very  sensitive  to  error.  Algebraic  errors  generally  lead, 
after  considerable  wasted  effort,  to  the  possibly  false  conclusion 
that  the  Lie  algebra  must  be  trivial. 

One  procedure  which  we  have  usefully  implemented  on  Macsyma 
adds  new  elements,  defined  equal  to  some  of  the  unknown  commutators. 

We  add  one  element  at  a  time,  substituting  all  newly  determined 
commutators  into  the  Jacobi  identities  as  we  proceed.  For  the  example 
which  we  will  describe  in  section  V  this  algorithm  led  to  a  closed 
Lie  algebra.  Furthermore,  when  this  algorithm  closes  the  Lie  algebra  it 
allows  us  to  determine  whether  a  free  parameter  such  as  that  associated 
with  the  eigenvalue  of  an  inverse  scattering  transform  exists. 
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For  those  sets  of  commutators  which  are  not  consistent  with 
any  nonabelian  closed  Lie  algebra,  the  algorithm  described  usually 
leads  quickly  to  an  Abelian  algebra.  An  example  of  this  is  the  WE 
analysis  for  the  combined  Burgers -KdV  equation 

u+  +  au  +  gu  +  vu  u  =  0.  (10) 

t  xxx  xx  '  x 


Taking 


4>x  =  P(u)$, 

the  WE  analysis  leads  systematically  to  a  set  of  commutation 
relations.  Applying  the  algorithm  described,  we  are  led  to  the 
conclusion  that  all  elements  of  the  Lie  algebra  must  commute.  This 
calculation  would  have  been  prohibitively  laborious  if  it  had  not 
been  done  on  Macsyma. 

The  algorithm  we  have  described  gives  a  useful  first  attempt  at 
closing  the  Lie  algebra.  For  the  example  we  present  in  section  V 
it  is  all  that  is  necessary.  But  the  algorithm  does  not  always 
terminate.  For  the  KdV  equation, 
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the  WB  ansatz 
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where 


[xrx2]  =  x3,  [Xj  ,x3]  =  x4>  [x2,x3]  =  x5 

[X2,X5]  =  [Xr5q  =0,  [X1,X4]  -  [X2,X]  (14) 

These  brackets  may  be  imbedded  in  arbitrarily  large  Lie  algebras. 

The  algebra  may,  however,  be  closed  without  adding  any  new  elements. 
Set 


x  =  c1x1  +  c2x2  +  ...  +  c5x5, 


(lb) 


where  the  CL's  are  constants  to  be  determined.  By  evaluating  all 
Jacobi  identities  we  determine  all  but  one  of  these  constants  and  all 


A 


10 


remaining  commutators.  One  constant  remains  arbitrary,  corresponding 
to  the  eigenvalue  in  the  associated  inverse  scattering  transform.  The 
ansatz  (15)  is  suggested  by  the  fact^  that  as  x+  ±»,  Q  T,  so  that 
X  must  contain  the  eigenvalue  if  we  are  to  solve  the  nonlinear  PDE  by 
an  inverse  scattering  method.  This  second  algorithm  appears  to  be  a 
very  general  one  for  finding  the  inverse  scattering  transform  of  an 
integrable  equation. 
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IV.  Implementation  on  Macsyma 

Macsyma  has  a  built-in  noncomnutative  product  operator  whose  purpose 
is  to  represent  multiplication  of  matrices.  For  our  purposes,  we  find  it 
convenient  to  use  this  noncomnutative  product  to  represent  the  Lie  bracket. 
Having  used  a  built-in  MACSYMA  operator,  we  can  apply  various  operations 
available  in  MACSYMA  for  evaluation  and  simplification  of  algebraic  expres¬ 
sions.  To  guarantee  that  the  Lie  bracket  will  be  handled  properly  by 
these  operations,  we  reset  several  of  the  MACSYMA  flags  associated  with 
the  noncomnutative  product.  By  doing  so  we  give  this  product  all  the 
desired  properties  except  antisymnetry,  which  must  be  built  into  the  table 
of  commutators. 

The  noncomnutative  product  operator  is  representated  in  MACSYMA  by  a 
dot,  and  so  is  called  the  "dot  operator."  The  double  bracket  which  arises 
in  the  Jacobi  identities 


[A,[B,C]]  , 

is  simply  represented  as 

A.(B.C)  . 

The  dot  operator  is  normally  taken  to  be  associative  in  MACSYMA.  By 
setting  the  DOTASSOC  flag  to  FALSE  we  make  it  nonassociative,  as  required 
by  the  Lie  brackets.  Similarly,  by  setting  the  DCTDISTRIB  flag  to  TRUE  we 
allow  A.(B+C)  to  simplify  to  A.B  +  A.C. 

To  distinguish  between  those  variables  which  are  elements  of  our  Lie 
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algabra  and  those  which  represent  numerical  coefficients,  MACSYMA  allows 
us  to  declare  variables  to  be  nonscalar.  This  is  done  using  the  function 
DECLARE.  The  way  in  which  scalars  are  handled  is  then  controlled  by  the 
flag  DOISCRULES.  We  set  DCTSCRUIES  to  TRUE.  If  c  is  a  scalar,  c.A  or 
A.c  will  then  simplify  to  c*A,  and  A.  (c*B)  will  simplify  to  c*(A.B) . 

With  DOTCONSTRULES  set  to  TRUE,  purely  numerical  coefficients  are  treated 
in  the  same  manner. 

The  table  of  commutators  is  stored  as  a  list  of  equations,  with  each 
equation  of  the  form 


A.B  •  C  .  (16) 

This  allows  us  to  use  the  EV  conmand,  which  will  evaluate  a  given  expression 
subject  to  a  list  of  equations.  For  example,  EV(A.B +D,TABLE)  yields  the 
expression  C+D  if  TABLE  is  a  list  of  equations  including  Eq.  (16).  If  EXPR 
is  a  more  complicated  expression  containing  Lie  brackets,  then  the  conmand 
EV(EXPR,TABLE)  will  substitute  C  for  A.B  each  time  the  latter  appears  in 
EXPR.  By  adding  the  argument  INFEVAL,  we  use  EV  in  an  "infinite  evaluation" 
mode.  If  D.C=E  Is  also  determined  by  TABIE,  and  D.(A.B)  appears  in  EXPR, 
EV(EXPR,TABI£,INFEVAL)  substitutes  an  E.  Note,  however,  that  EV  does  not 
distribute  over  the  dot  product.  If  A.B  =E+F,  then  only  the  first  layer  of 
D.  (A.B)  can  be  evaluated  by  an  EV  conmand.  To  evaluate  further  we  must 
first  use  the  EXPAND  conmand.  This  conmand  yields  a  simplified  expression 
in  which  the  dot  and  scalar  products  have  been  distributed  over  sums  of 
terms,  scalar  multiplication  has  been  factored  out  of  the  dot  product,  etc. 


L  k 


A 


The  program  reads  in  the  number  of  initial  elements  (N)  in 
the  unclosed  algebra,  and  then  reads  in  the  known  structure  constants 
of  the  Lie  algebra,  using  the  READ  conmand  in  a  loop.  Of  course  at 
most  N(N-l)/2  such  constants  must  be  given,  with  the  antisymmetric 
twin  of  each  equation  automatically  added  to  the  table,  and  with  A.  A 
automatically  set  to  zero.  Given  the  structure  constants,  the  program 
proceeds  to  evaluate  all  Jacobi  identities,  using  two  layers  of 
EXPAND  and  EV  commands  for  each. 

To  implement  the  algorithms  described  in  the  previous  section, 
the  program  defines  functions  which:  evaluate  a  given  Jacobi  identity; 
update  all  Jacobi  identities  using  new  information;  update  a  given 
Jacobi  identity;  add  a  given  equation  to  the  table;  evaluate  all 
Jacobi  identities  involving  a  newly  added  element  of  the  algebra;  store 
the  updated  table  along  with  all  defined  functions,  options,  and 
properties  on  disk;  display  the  table  of  structure  constants.  These 
functions  perform  the  laborious  algebraic  tasks  involved  in  trying 
to  close  the  Lie  algebra.  They  allow  the  user  to  implement  the  algo¬ 
rithms  of  the  previous  section  in  a  flexible  fashion,  making  use  of  any 
new  insight  as  he  gains  experience  in  attempting  to  close  such  an  algebra. 
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V.  An  Example:  The  Kinetic  Three-Wave  Equations 

We  have  used  Macsyma  to  close  the  Lie  algebra  for  Eqs.  (1). 

Here  we  first  derive  the  appropriate  commutation  relations,  and  then 

describe  our  use  of  Macsyma  to  close  the  algebra. 

Eqs.  (3)  are  our  starting  point,  with  P  and  Q  each  assumed  to 

depend  only  on  the  three  N. 's.  Setting  the  coefficient  of  N.  equal 

J  Jx 

to  zero  for  each  j ,  we  obtain 

VjPN'  •  ,  j  -  1,2,3.  (17) 

Taking  cross  derivatives  of  these  equations,  we  find 


(18) 


Since  the  group  velocities  are  assumed  all  different,  Eq.  (18)  implies 


(19) 


or 


3 

P-  l  g.PU+Z,  (20a) 

j-1  3  3 

3 

Q  *  -  l  v.g.(N.)  ♦  W.  (20b) 

j-1  J  J  J 


0,  so  that  gj  must  be 


From  Eq.  (7)  we  further  find  that  g V '  = 
of  the  form 


g.  =  N.X. 
S3  J  3 


y 


where  we  have  absorbed  the  possible  constant  teim  into  W  and  Z. 
The  Jacobi  identities  imply  that  each  must  commute  with  all 
elements  of  the  algebra.  So  the  VVs  can  be  set  to  zero  without 
changing  the  structure  of  the  algebra. 

We  have  concluded  that  P  and  Q  must  be  of  the  form 
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P  =  l  X  N  ♦  Z, 
j»l  3  J 

(21a) 

3 

q  3  -  y  V,X.N  +  W. 
j»l  3  3  J 

(21b) 

Substitute  these  expressions  back  into  Eq.  (7),  setting  the  coefficient 
of  each  term  in  the  resulting  polynomial  in  the  N. 's  to  zero,  to 
obtain  a  set  of  commutation  relations.  To  simplify  the  expression 
of  these  relations  we  define 

\  -  Xr  Y2  =  -  X2,  Y3  =  X3.  (22) 


In  terms  of  the  Y^’s,  our  commutation  relations  are 


1 6 


[Y.Wj  s  v JZ.YJ, 
m  m  m 


[W,Z]  =  0. 


The  commutator  of  Y„  with  Z  or  W  is  not  determined, 
m 

Now  we  apply  the  first  algorithm  described  in  section  III, 
using  Macsyma.  Let 


[Z,Yj]  =  Ur 


Evaluating  all  Jacobi  identities,  we  find  that  the  commutators  of  Uj 
with  W  and  Z  are  now  left  undetermined.  Let 


[Z,U.]  =  0.. 


Evaluating  all  Jacobi  identities  we  find  that  we  must  again  extend 
the  algebra.  Let 


[Z,U.]  =  U. 


Now  when  we  evaluate  the  Jacobi  identities  the  algebra  collapses  back 
down.  We  find  that  IL  =  0  for  all  j.  We  obtain  relations  between 
various  of  the  elements.  Substituting  back  into  the  Jacobi  identities, 
we  are  able  to  determine  all  commutators . 
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We  reduce  the  number  of  elements  in  our  Lie  algebra  by  setting 
to  zero  all  elements  or  combinations  of  elements  which  commute  with 
every  other  element.  This  does  not  affect  the  conmutation  relations, 
but  does  simplify  our  search  for  a  linear  representation. 

As  a  final  simplification,  we  work  in  a  reference  frame  where 
some  of  the  structure  constants  vanish.  We  choose  our  reference 
frame  so  that 


vl|v3  -  v2l  +  v2|v1  -  v3l  ♦  v3|v2  -  vx|  =  0.  (24) 


Since  we  have  assumed  that  all  three  group  velocities  are  different, 
we  can  always  find  such  a  frame.  In  this  reference  frame  we  can 
close  the  Lie  algebra  by  setting 


sgnCv7  -  v  )  +  sgn(v,  -  v.) 
[Y1’Z1 - VX(V3  -  -  v/  •  W  • 


(25a) 


sgn(v,  -  v7)  +  sgn(v,  -  v.) 

[y5'z) - M  - 


(25b) 


Y 


2 


<V1 


-  V,)  (: 


3^vl 


(25c) 


These  relations,  along  with  Eqs.  (23),  define  a  closed  Lie  algebra. 

The  final  result  looks  deceptively  simple.  It  probably  would 
have  taken  a  good  deal  of  guessing  and  playing  around  with  the 
algebra  to  arrive  at  Eqs.  (24)  and  (25)  without  the  systematic  approach 
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we  used.  And  the  systematic  approach  would  have  been  discouragingly 
laborious  without  the  use  of  Macsyma.  Furthermore,  having  gone 
through  our  systematic  procedure,  we  are  now  assured  that  the  algebra 
we  obtained  is  the  largest  possible.  This  is  of  interest  because  it 
implies  that  there  can  be  no  free  parameters  of  the  sort  that  arises 
when  a  nonlinear  PDE  can  be  solved  by  an  inverse  scattering  method. 

Using  the  closed  Lie  algebra  thus  obtained,  we  have  found  a 
linear  3x3  representation,  and  have  thus  been  able  to  imbed  the 
nonlinear  Eqs.  (1)  in  the  linear  Eqs.  (3).  We  have  been  able  to  use 
this  imbedding  to  obtain  an  interesting  set  of  exact  multi-shock 
solutions.  Except  for  the  single  shock  solutions  these  solutions  are 

7 

new.  They  are  described  elsewhere  . 
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